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The role of Supernova Remnants (SNRs) for the production of 
the Galactic Cosmic Rays is reviewed from the point of view 
of theory and very high energy gamma-ray experiments. The 
point is made that theory can describe young SNRs very well, 
if the evidence from the synchrotron emission is used to em- 
pirically determine several parameters of the theory, and thus 
theory can predict the relative contributions of hadronic and 
leptonic gamma rays at TeV energies. This is exemplified for 
several objects that have been observed intensively during the 
last years. Future key observations are discussed. 



1 Introduction 

Shell-type supernova remnants (SNRs) are widely assumed to be 
the sources of the Cosmic Rays (CRs), as they are observed in the 
neighborhood of the Solar System. This concerns particle energies 
up to the "knee" of the energy spectrum at several 10^^ eV or 
possibly beyond (see P for a recent review). From estimates of 
the Galactic Supernova (SN) rate and the CR escape rate from 
the Galaxy SNRs have then to convert on average about 10% of 
their entire mechanical explosion energy into CRs - an enormous 
requirement. 

A direct experimental investigation of SNRs as CR sources is 
possible with 7-ray observations at very high energies > 100 GeV 
(VHE). The argument is the following: the acceleration of parti- 
cles to CR energies is assumed to occur primarily at the outer, 
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Figure 1: ASCA image of SN 1006 
in hard X-rays from |2]. The emis- 
sion comes mainly from the North- 
eastern and Southwestern areas, in- 
terpreted as polar caps, where the 
mean magnetic field is quasi-parallel 
to the shock normal of the outer SNR 
blast wave The resolution of 

ASCA is slightly better than that of 
the H.E.S.S. array, indicating what is 
presently achievable in VHE 7-rays. 



Figure 2: Schematic of the average 
magnetic field geometry for a SN 
explosion into a uniform interstellar 
medium with a homogenenous field 
—1- ®nB angle between the 

shock normal vector n and B^i . Injec- 
tion of downstream suprathcrmal ion 
is only possible for sufficiently small 
values of i.e. in the hatched 

polar regions. 



quasi-spherical shock which also compresses and heats the ambi- 
ent circumstellar medium. The expanding shock wave confines the 
accelerated particles in its interior until its velocity decreases sub- 
stantially at late times. Then the shock gets "old" and the more 
energetic particles successively leave the remnant. There fore a dis- 
tant 7-ray observer of a "young" SNRs can see the 7r°-decay and 
nonthermal Bremsstrahlung (NB) emission, due to CR collisions 
with gas particles in the interior, jointly with the Inverse Comp- 
ton (IC) radiation as originating from a localized source. Together 
with the electron synchrotron spectrum - from radio to hard X-ray 
energies - and the synchrotron morphology (e.g. Fig. P) this is 
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Figure 3: Cas A in X-ray syn- 
chrotron light, observed with 
Chandra. Image courtesy of 
NASA/CXG/GSFC/U. Hwang et 
al.. 
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Figure 4: Synchrotron spectral en- 
ergy density (SED) of Cas A: data 
and theoretical spectrum The 
low-frequency photon spectrum has a 
spectral index a « 0.8 (dotted curve), 
not quite increasing to the test parti- 
cle value a ~ 0.5 before in the sub- 
mm region the electrons inside the 
SNR are already radiatively cooled, 
so that a K, \ before the cutoff. 



the nonthermal electromagnetic evidence. I want to argue below 
that the synchrotron emission allows us to seperate the contribu- 
tions of nuclear CRs and ultrarelativistic CR electrons to the 7-ray 
emission; in addition, we can then determine key theoretical quan- 
tities like the effective magnetic field and the rate of injection of 
suprathermal ions into the acceleration process. In this way it be- 
comes possible to calculate the energy density and the energy spec- 
trum of nuclear CRs in the SNR from theory. To obtain quantita- 
tive results, obviously some further astronomical multi-wave-length 
(MWL) information is required. This concerns the source distance, 
the angular size/morphology, the SNR expansion rate, and the cir- 
cumstellar density structure. SNe also result from two basically 
different physical processes whose main manifestations are type la 
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explosions (deflagration/detonation of an accreting White Dwarf) 
and type II explosions (core collapse of a massive star), with sev- 
eral variants according to the mass of the progenitor star. These 
different SN types are connected with different magnitudes of the 
ejected mass. Therefore also the explosion type must be known 
from astronomical measurements. 

In this talk I will review the general status of the 7-ray obser- 
vations, emphasizing the non-spherical aspects of the nonthermal 
emission from SNRs and the relation of theory and experiment in 
this field. In this light I will describe four major sources: Cas A, 
SN 1006, RX J1713. 7-3946, and Vela Jr. I will end with a short 
discussion of the contribution to nonthermal SNR research by the 
recent Galactic Plane Scan, performed with the H.E.S.S. experi- 
ment. In some respects this is an extension of a paper given at the 
28th ICRC in Tsukuba j4j. 

2 Gamma-ray detectability of SNRs and 
their non-spherical aspects 

The EGRET instrument on the Compton Gamma Ray observatory 
(GGRO) has not been able to find unequivocal evidence for 7-ray 
emission from SNRs in its energy range below a few GeV. The 
reasons are the low source flux, even for objects as close as 1 kpc, 
and the large angular extent of ~ 1° - of the order of the diameter 
of the full Moon - which implies a high 7-ray background from 
diffuse Galactic CRs. If we anticipate the particle energy spec- 
tra in the sources to be much harder than the energy spectra of 
the diffuse Galactic CRs, then the signal-to-background ratio de- 
creases with decreasing 7-ray energy, making the GeV range gener- 
ally unfavorable for detection. In contrast, ground-based imaging 
Cherenkov telescopes have made several detections in dedicated 
VHE observations (see section 4), and the H.E.S.S. experiment 
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has identified SNR counterparts for several sources found in its 
Galactic Plane Survey. The reasons are the large effective area 
of these telescopes/telescope systems, the relatively lower diffuse 
7-ray background at TeV energies, and the much better angular 
resolution of ~ 10~^ degrees compared to EGRET. This makes 
nearby SNRs at least marginally detectable 

2.1 Non-spherical aspects of SNRs 

Theoretical models for diffusive shock acceleration at SNRs face the 
difficulty of having to cope with the fundamentally non-planar and 
even non-spherical geometry of a point explosion into an environ- 
ment that lacks spherical symmetry. The dynamics is described by 
kinetic equations for the particle distributions f{p, r, t) as functions 
of particle momentum p, radial distance r and time t, nonlinearly 
coupled with the hydrodynamics of the thermal gas. Only spheri- 
cally symmetric solutions are available until now which solve this 
intrinsically time-dependent problem [HIIZI- It is clear on the other 
hand that the magnetic field, which regulates the particle injection 
rate into the acceleration process, cannot be spherically symmetric 
and is even on average at best axially symmetric in SNRs. 

The simplest case is a type la SN in a uniform interstellar 
medium and magnetic field, with SN 1006 as the clearest exam- 
ple. (Fig. P). The time-average magnetic field line geometry is 
schematically shown in Fig. |21 

For kinematic reasons injection of suprathermal ions escaping 
from a thermalized downstream region can only occur for quasi- 
parallel shocks, where the instantaneous angle ©j^g ^ 7r/2. And 
clearly acceleration can occur only at those parts of the shock sur- 
face, where particles can be injected. Particle acceleration is also 
directly connected with the self-excitation of Alfvyen waves which 
stochastically change ©^g. As a consequence we have (i) a stochas- 
tic self-limitation of the ion injection rate rj through nonlinear wave 
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production, from rj\\ ^ 10^^ to an T^gg 10^^, plus (ii) a system- 
atic reduction of rj due to the overall average field morphology, i.e. 
strong wave production can occur only locally in the polar regions, 
(iii) the hadronic 7-ray emission is therefore also dipolar for uni- 
form external field B_i, and (iv) the same is true for the synchrotron 
emission as a result of field amplification by factors between 5 and 
10 in the ion acceleration regions ^2113^1! with essentially lower 
emission from the extensive equatorial region. This last conse- 
quence has been impressively proven in a recent analysis of the 
XMM data for SN 1006 by [TT]. 

Altogether this injection asymetry requires a reduction of the 
overall acceleration efficiency of nuclear particles as calculated in 
the spherically symmetric model. The reduction factor is given by 
the ratio of the polar areas to the total shock surface area. This 
ratio is about 0.2 for like SN 1006. In order to reach an 

overall acceleration efficiency of 10% this requires the shock regions 
in which acceleration actually occurs to achieve an acceleration 
efficiency of about 50%. Such a high efficiency implies an extremely 
nonlinear acceleration process with a strong backreaction of the 
accelerated particles on the shock structure. 

3 Comparison with theory 

The comparison with theory is of course an essential aspect. How- 
ever, at present the theory is still incomplete. The full solution 
of the Fokker-Planck transport equations for the distribution func- 
tions f{p,r,t) of nuclear particles and electrons, coupled with the 
hydrodynamics of the thermal plasma through the CR pressure gra- 
dient and wave dissipation, requires even in spherical symmetry the 
knowledge of several "unknowns" : the effective, amplified magnetic 
field strength -Bgff, the actual proton injection rate ?7gg, and the 
amplitude of the electron distribution. These unknowns can only 
be determined through an analysis of the synchrotron observations 
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Figure 5: Data of an individual 
Chandra 2-10 keV filamentary struc- 
ture in Cas A and model fit to these 
data, interpreted as the result of 
strong postshock synchrotron losses 
(at ?/' < 0). 



Figure 6: Gamma-ray spectral en- 
ergy density for Cas A. The (uncon- 
firmed) HEGRA detection and the 
EGRET upper limit are shown to- 
gether with the theoretical predic- 
tion for the tt'^ -decay (full line), IC 
(dashed line) and NB (dash-dotted 
line) emissions. 



which involve the relativistic electron component. This analysis is 
actually possible, because for particle energies E ^ mpc^, corre- 
sponding to ultra-relativistic protons, electrons behave like protons 
in the acceleration process (e.g. [T^l. 

I will only summarize the situation here. 

The electrons are parasitically accelerated in an environment 
produced by the accelerating nuclear particles because they can 
not modify the shock themselves due to their small energy density. 
And therefore at energies E ^ mpc^ their momentum distribution 
equals in form that of the nuclear particles except for radiative 
(synchrotron) losses: 

• The radio synchrotron spectrum is generally steeper than in 
the test particle approximation, because the radiating low- 
energy electrons "see" only the discontinuous subshock in the 
thermal gas and not the full shock transition that includes 
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the extended CR precursor. Interpreting this effect in terms 
of nonhnear shock modification by the accelerated nuclear 
component determines the ion injection rate rj^Q 

• The energy of the radio electrons should be substantially 
lower than lOmpc^. Together with the requirement to fit the 
entire synchrotron spectrum, including the cutoff at hard X- 
rays, this determines the (amplified) magnetic field strength 

^eff 

• The amplitude of the relativistic electron density then follows 
from the amplitude of the synchrotron spectrum, which fixes 
the electron-to-proton ratio in the accelerated CRs 

4 Individual SNRs 
4.1 Cas A 

Cas A is presumably the result of a so-called Type lb SN, the 
core collapse of a massive Wolf-Rayet star that has already shed 
its hydrogen envelope through a fast stellar wind (Fig. Oj). In a 
detailed model for the thermal X-ray emission by (14j the final 
Wolf-Rayet wind phase has compressed the inner part of the dense 
slow wind from the preceeding Red Supergiant (RSG) phase into a 
dense luminous shell; in turn the subsequent SNR shock has already 
reached the unperturbed RSG wind region beyond the shell. 

The strong deviation of the synchrotron spectral energy density 
from a test particle spectrum below some tens of GHz (Fig. E]) 
implies a strong modification of the shock by accelerated nuclear 
particles, an amplified post-shock field B^q ^ 200yuG, and a field 
Sgg ^ 500/iG in the RSG wind shell dSJ. 

The X-ray morphology of Cas A from Chandra also shows strongly 
pronounced filamentary structures of the outer shock (Fig. E)), an- 
alyzed by [in] and [T7j . The multi-Te V electrons accelerated at the 
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shock form a very thin quasi-spherical shell - thinner than that of 
a children's rubber ball - as the result of violent synchrotron cool- 
ing. This cooling scale determines the interior amplified field. It 
turns out that this field "measurement" agrees with that using the 
spectral distortion in the radio frequency range within the errors of 
20 to 30 percent. Similar results have been obtained for SN 1006 
da CHI and Tycho's SNR 

The theoretical prediction of the 7-ray fluxes (Fig.lHl from jl5j ) 
shows a TT^-decay 7-ray flux that dominates those from IC scat- 
tering and NB by two orders of magnitude at 1 TeV, making it a 
clear hadronic 7-ray source by a large margin. This prediction es- 
sentially agrees with the flux detected by the HEGRA experiment 

EH- 

The IC and NB flux determinations are quite robust results. 
A reliable independent measurement of the 7-ray flux with the 
large Northern Hemisphere telescopes VERITAS or MAGIC would 
therefore be of paramount importance. 

4.2 SN 1006 

This SNR has been observed in hard X-rays with ASCA to show 
purely nonthermal emission from the two hot spots at the poles, as 
discussed before, and this emission was interpreted as synchrotron 
radiation j2]. Later high- resolution C/ianc?ra and XMM observations 
strengthened this picture. However, the TeV 7-ray detections by 
the single CANGAROO telescopes CANGAROO I ^ and CAN- 
GAROO II [231 could not be confirmed in a total of 24.5 hours of 
observation time by the H.E.S.S. stereoscopic system |21]. Recent 
CANGAROO stereo observations could not detect the source any 
more either and have led to the withdrawal of the earlier 7-ray 
detection claims P^. 

There are two reasons for the 7-ray non-detection. First of all, 
the magnetic field in the SNR interior is considerably amplified 
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Figure 7: Integral VHE photon flux 
from the northeastern polar cap of 
SN 1006. The value Bq = 30^ G cor- 
responds to an interior effective field 
strength of ISO^G. The CANGA- 
ROO and H.E.S.S. data are shown 
together with the theoretical flux esti- 
mates in ]iy!/ . 



Figure 8: H.E.S.S. image of RX 
J1713.7-3946 at-f-ray energies S 800 
GeV, with the ASCA 1 - 3 keV X-ray 
contours superposed 120^ . 



(Sgg ~ 150/iG), SO that the IC radiation is strongly suppressed, 
given the observed synchrotron emission. Secondly, the external 
hydrogen density A'^^^ is in all probability quite low, A'^^^ < O.lcm"^. 
Since in the Sedov phase, in which SN 1006 is at present, the vr"- 
decay 7-ray flux is proportional to N]j, the low gas density 
implies a low hadronic 7-ray emission as well. This situation has 
been analyzed in detail by JH]- Given the lowest value of Nh = 
0.05 cm~^, discussed in the literature, we expect the 7-ray emission 
to be only a factor of 3 smaller than the present H.E.S.S. upper 
limit for the northeastern polar cap (Fig. [7j). Since SN 1006 is the 
simplest case of a SNR with strong nonthermal emission in radio 
and X-rays and therefore the theoretically best understood object 
of its kind, it would obviously be important to detect it also in TeV 
7-rays in a deep observation of about 200 hours with H.E.S.S. or 
CANGAROO. 
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4.3 SNR RX J1713.7-3946 



Originally found in the ROSATX-iiiY survey [26j, SNR RX J1713.7- 
3946 was detected at VHE energies by CANGAROO ^ and inter- 
preted as an IC-dominated TeV source. ASCA observations PH} !^ 
had shown that the X-ray spectrum is a purely nonthermal con- 
tinuum. In 2002 the CANGAROO group revised its interpretation 
and rather favored a hadronic scenario from the shape of the TeV 
energy spectrum jSO]. This phenmenological interpretation was 
questioned by [HI] and jH2] on the grounds that the upper limit 
from a nearby EGRET source 3EG 1714-3857, believed to be as- 
sociated with the TeV source, was inconsistent with the hadronic 
extrapolation of the GANGAROO spectrum to the GeV region. 
Such different views were not entirely surprising at the time, con- 
sidering the complex morphology of the general region in which 
the source is embeded. This is particularly well visible in the CO 
map. Nevertheless, at the present time CO data jHSl give a most 
probable kinematic source distance of 1 kpc. 

The H.E.S.S. experiment subsequently confirmed the GANGA- 
ROO detection and could give for the first time a spatially resolved 
VHE image of a SNR [33], whose overall shell structure correlated 
closely with the shell structure in ASGA 1-3 keV X-rays. The 
remnant diameter is about 1°. This was unambiguous proof for 
the acceleration of charged particles to energies beyond 100 TeV^ 

Apart from the fit shown in Fig.El the H.E.S.S. VHE differential 
spectrum can be equally well fitted by a power law with exponential 
cutoff oc E-^exp--^/^^ with T = 1.98 ± 0.05 and = 12 ± 2. 
At energies E <^ Ec this measured spectrum corresponds to the 
test particle limit of diffusive shock acceleration theory for nuclear 
particles in a strong shock. Even ignoring nonlinear backreaction 
effects the extrapolation of the charged particle spectrum to lower 

^The following discussion of this source and the figures |H1 and El are based on 
the more recent H.E.S.S. results, not yet released at the time of the conference. 
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energies with a proton spectrum oc E^'^ gives a hadronic 7-ray 
spectrum below an improved upper limit of EGRET which now 
assumes that RX J1713. 7-3946 is not linked to the known EGRET 
source 3EG 1714-3875 |35]. 

Making in fact the best case for an IC interpretation by assum- 
ing a very low magnetic field strength of about lO/xG inside the 
SNR, the resulting IC spectrum fits the H.E.S.S. data quite poorly. 
This outweighs the good correlation between the X-ray synchrotron 
and the 7-ray emissions which at first sight would suggest a lep- 
tonic origin of the 7-ray emission as well. In addition, already a 
small amplification of the magnetic field in the remnant rules out 
a dominant IC emission and a fortiori a dominant NB, whereas a 
hadronic 7-ray spectrum continues to fit the data quite well jSSI. 
Despite the complex CO morphology and unknown age of this rem- 
nant I believe that a hadronic interpretation of the 7-ray emission 
from RX J1713. 7-3946 is clearly favored. The most plausible accel- 
eration scenario is that of a massive progenitor star which, over its 
long evolution time, had produced a large stellar wind bubble into 
which finally the SN exploded about a thousand years ago (see e.g. 
Fig. 5 of [SEl)- Only this makes the almost circular X-ray morphol- 
ogy understandable that was found by ASCA and XMM [211 IHHj . 
And it can explain the low effective density inside the SNR |^ . A 
detailed theoretical model is needed to investigate the consistency 
of such a picture with the existing MWL evidence. 

4.4 Vela Jr. 

The SNR RX J0852.0-4622, also called Vela Jr., was also found with 
ROSAT as a very large 2° x 2° quasi-circular X-ray structure 
and later confirmed by ASGA [501 In the VHE range it was 
detected by GANGAROO g2] and by H.E.S.S. gg, see Fig. CHI 

The flux above 1 TeV equals 1.4 times the flux from RX J1713.7- 
3946 and is about equal to the flux from the Crab Nebula. The 
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Figure 9: H.E.S.S. differential 7- 
ray spectrum of RX Jl 713. 7-3946 for 
the whole region of the SNR solid 
black circles. The best fit of a power 
law luith energy- dependent photon in- 
dex is plotted as a black line. The 
H.E.S.S. 2003 data are given by the 
blue open circles. Error bars are ilo" 
statistical errors fSSj/ . 
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Figure 10: Vela Jr. in TeV j-rays. 
The color scale is in number of ^ -ray 
events. The total significance of this 
3.2 hours lifetime observation with 
H.E.S.S. is 12(7. The source radius 
is four times the radius of the full 
Moon. 



differential spectrum is a hard power law spectrum with F = 2.1 ± 

0.1 ±0.1 gn]. 

The physical characteristics of Vela Jr. appear to be rather 
similar to those of RX J1713. 7-3946, even though the source may be 
even closer and younger. Interestingly, also a narrow filament was 
recently found in Chandra data |13] . If we interpret this filament 
as part of the outer SNR shock its sharpness suggests magnetic 
fields at least in the lOOyuG range, almost independently of the 
detailed astronomical properties of the source [IHI- Therefore again 
a hadronic interpretation of the 7-ray emission is plausible, even 
though little theoretical analysis has been performed up to now. 
Further and more detailed results on the morphology and spectrum 
are expected to be available soon. They should make this source 
the second major VHE SNR source in the Southern Hemisphere. 
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5 Galactic Plane Scan SNRs, Conclu- 
sions 



The H.E.S.S. Galactic Plane Scan |lBlll7j has revealed some 20 odd 
new sources until now. Several of these soures can be associated 
with SNRs with reasonable certainty, but many of them are still un- 
identified at the present time. It is clear that further study of these 
SNRs is required, both at VHE energies and at other wavelengths. 
An example is HESS J1813-178 where, quickly following the initial 
H.E.S.S. publication of the first scan results, a SNR was found in 
existing VLA radio data f^, and a known coincident ASCA source 
was found also by INTEGRAL in the 20-100 keV band |49j. 

In the present context such studies are important to identify the 
statistical VHE properties of the Galactic SNR and thus to control 
whether the few very nearby objects, that can be studied in great 
detail, are typical representatives of this population or not. 

Even though it is difficult to make such population studies in 
the Northern Hemisphere, there are at least two bright historical 
SNRs in this Hemisphere, Gas A and Tycho's SNR. Gas A has 
been detected by HEGRA in TeV 7-rays with rather low signifi- 
cance, whereas the HEGRA fiux upper limit for Tycho's SNR jSOj is 
only a factor of unity above the theoretical VHE prediction jSUIHl- 
Together with other known Northern Hemisphere SNRs, the con- 
firmation/detection - or not - of these objects with VERITAS or 
MA GIG would be of high significance for the question of GR origin 
in SNRs. 
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